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O ■ Abstract 

(N ■ 

We analytically compute the dominant contributions to the /3-functions for the top- Yukawa 
coupling, the strong coupling and the Higgs self-coupling as well as the anomalous dimensions 
of the scalar, gluon and quark fields in the unbroken phase of the Standard Model at three- 
■ loop level. These are mainly the QCD and top- Yukawa corrections. The contributions from 

the Higgs self-interaction which are negligible for the running of the top- Yukawa and the 
r S^ ' strong coupling but important for the running of the Higgs self-coupling are also evaluated. 

Q 1 Introduction 

£SJ ' Using perturbation theory in any renormalizable quantum field theory comes with the price 

that the parameters of this theory, e.g. the couplings and masses, in general depend on the 
, renormalization scale //. The precise description of the evolution of these parameters with the 

energy scale is an important task in any model. This is done by means of the Renormalization 
Group functions, that is the /3-functions and anomalous dimensions. The knowledge of the 
three-loop contributions to the /3-functions for the Standard Model (SM) and its extensions 
is important for physics at the very high energy frontier and for cosmology. Here are some 
examples: The running of the gauge couplings plays an important role for the construction 
of Grand Unified theories of the strong and electroweak interactions. The /3-functions for 
the scalar self-interaction and for the top quark Yukawa coupling constant are important for 
the analysis of Higgs- inflation in the SM [IH1]- The current investigations of these issues are 
based on the two-loop approximation. The inclusion of the next order could be essential. 
Further, in a recent work [3], the possibility has been discussed that the SM, supplemented 
by the asymptotically safe gravity could play the role of a fundamental, rather than effective 
field theory. Within this framework the mass of the Higgs boson has been predicted to be 
approximately 126 GeV. The theoretical uncertainty of the prediction is about 2 GeV. 

Recent exciting evidence from several SM-like Higgs search channels at both the CERN 
Large Hadron Collider and the Fermilab Tevatron [BHS] point to the possibility of a SM 
Higgs boson with a mass in the vicinity of 125 GeV which is in truly remarkable agreement 
with the aforementioned prediction^] This calls for more precise calculations, in particular of 
/3-functions, in the SM. In the present paper we are particularly interested in the evolution 
of the Higgs self-coupling as well as the top- Yukawa coupling in the SM. 



Note that the boundary condition \(Mpi an ck) = 0, leading to the prediction of the Higgs mass close the 
the experimental evidence, has been also discussed recently in [5]. 
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The underlying gauge group of the SM is an SU C (3) x SU(2) x U y (l) which is spontaneously 
broken to SU C (3) x Uq(1) at the electroweak scale. As the renormalization constants for 
fields and vertices do not depend on masses and external momenta in the MS-scheme, we will 
perform our calculations in the unbroken phase of the SM. 

The most important contributions to the running of the Higgs self-coupling A arise from 
the top- Yukawa coupling and the strong sector. All other Yukawa couplings are signif- 
icantly smaller due to to the smallness of the respective quark masses. From the top 
mass M t ~ 172.9 GeV we get the Yukawa coupling at this scale y t (M t ) = \/2^ « 1 where 
v ~ 246.2 GeV is proportional to the vacuum expectation value of the scalar field $ from 
which results the Higgs field after the spontaneous symmetry breaking: |($)| = ^j=. 

The next Yukawa coupling to be considered would be y& = \f2^ ps 0.02. The strong cou- 
pling at the scale of the Z boson mass is g s {Mz) ~ 1-22 whereas the electroweak couplings 
gi(Mz) = cotdw ~ 0-36 and gi{Mz) = ~ 0.65 give much smaller contributions which 

are further suppressed by the isospin and hyper charge factors. For this reason we will con- 
sider a simplified version of the SM or - as one could also see it - a minimal extension of QCD 
by setting <?i = <?2 = in our calculation. For a Higgs mass of 125 GeV the value of the Higgs 
self-interaction would be X(Mh) ~ 0.13 at the scale of the Higgs mass. The relevance of this 
parameter will be examined in section [21 

The outline of the work is as follows. In the next section we discuss the main definitions 
and the general setup of our work. Section deals with the technical details, including the 
treatment of 75 . In sections 0] and we present our results for the /3- functions of the top- 
Yukawa, the strong and the Higgs self-couplings and the relevant field anomalous dimensions. 
The numerical influence of the computed three-loop corrections on the evolution of the quartic 
Higgs coupling is discussed in section For this analysis we will include the already known 
contributions with g\ and gi at one-loop and two-loop level. Finally, section [7] contains our 
conclusions and acknowledgements. 

All our results for /3-functions and anomalous dimensions can be retrieved from 
http: / /www-ttp. particle.uni-karlsruhe.de /Progdata/ttpl2/ttpl2-012/ 



2 General Setup 

The Lagrangian of our model consists of three pieces: 

£ = C-QCD 4" C vt + (1) 

The QCD part is defined by 

Cqco = - \G%G a ^ - {d^A^f + d^c a d»c a + gj abc d^A^t? 



q ^ 



(2) 



where q runs over all quark flavours, the gluon field strength tensor is given by 

G% = d^Al - d v Al + gJ abc A^A c v (3) 
and f abc are the structure constants of the colour gauge group with the generators T a : 
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The complex scalar field <1> and the left-handed parts of the top and bottom quarks t L and 
b L are doublets under SU(2): 

-(£)■ *-(:).■ 

Setting all Yukawa couplings to zero except for the top coupling y t the Lagrangian for the 
Yukawa sector is given by 



C vt = -y t {i R & c Q L + Q L $ c t R } 

= -y^2,-*i)-(l) L +^ L -(%)t R } (6) 

= -y t {(tP R t) <S>* 2 + (tP L t) $ 2 - (bP R t) - (fP L 6) $1} . 
Finally, we have the scalar sector of the model 

= d^d^<S>-m 2 &<S>- A(V$)\ (7) 



The indices L and R indicate the left- and right-handed part of the fields as obtained by the 
projectors 

Pl = \{1-1s) P R = l -{l + i,). (8) 



This model is renormalized with the counterterm Lagrangians 

5C QCD = - l -5zf 9) (d,A a v - d u Alf - Uzf 9) g s f hc (d,A a u - d u A«) A^A 



V 



- \bzf 9) g 2 a {j ahc A\A c ^ + SZ^d^d^c 11 + 5z[ CC9) g s f abc d^c a A b ^c c 
+ E {{ft [ 5Z t L ] P L + 6Z$P R ] q + 9 -qA a T a [SZ[ 999) P L + Sz[ 9 fP R 



(9) 



q 

for the QCD part, 

SC Yukawa = -5z[ m y t { (tP R t) $* 2 + (iP L t) <S> 2 - (bP R t) $J - (tP L b) $1} (10) 

for the Yukawa part and 

5£* = SZ^d^d^ - m 2 5Z^&$ + 5Z[^ ] (Vd>) 2 (11) 

for the ^-sector. Note that in general the left- and right-handed parts of quark fields and 
quark-vertices are renormalized differently!! So the renormalization constant for the strong 
gauge coupling g s can be obtained for example from 

r?(ttg) ^(ttg) 

^7 b^ !i? (12) 

93 ~ 7 (2t) rJlgj ~ 7 (2t) nj2g~) 



or the renormalization constant for y t from 



z {tb<£) 



Z Vt = 1 . (13) 



7 (2t) y (2t) 7 (2<K) 
^2,L ^2,Ft. ^2 



2 In our case this is true for the quark fields participating in the Yukawa sector. As we do not consider the 
electroweak interaction here and neglect all Yukawa couplings except for y t the light quark fields u,d,s and c 
have the same renormalization constant for the left- and right-handed part. 
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Here the renormahzation constants have been defined in a minimal way as 



Z = 1 + 5Z, (14) 

with 5Z containing only poles in the regulating parameter e = (4 — D)/2 of the dimensional 
regularization and D being the engineering space-time dimension. The Higgs self-coupling A 
is related to the Higgs mass at tree level Mh via 

which for M H = 125 GeV yields X(M H = 125GeV) « 0.13. For the running of the top- 
Yukawa coupling the contribution from A is negligible compared to the top- Yukawa and strong 
coupling. The corresponding four-$ vertex is nevertheless needed for the renormahzation at 
three-loop level, namely to kill the subdivergence from the fermion loop in diagrams like Fig. 
Q] (d) . The j3- function for a coupling X is defined as 

and is given as a power series in all couplings of the model, namely g a , y t and A. Note that 
the /3-functions (5 gs and f3 Vt are proportional to g s and y t respectively whereas f3 x has one 
part proportional to A and one part proportional to yf with no A-dependence at all. The 
anomalous dimension of a field / is defined as 

- dkiZj 1 _™ 1 /(n) 

where Z f is the field strength renormahzation constant for the respective fieldH The j3- 
functions for all couplings are independent of the gauge parameter £ whereas the anomalous 
dimensions of the fields are not. 



3 Calculation 

As we are interested in /3-functions it is enough to compute the UV divergent part of all 
diagrams in order to determine the necessary renormahzation constants. In the MS-scheme 
the latter depends only polynomially on external momenta and masses. Therefore most of our 
renormahzation constants could be computed from massless propagator-like diagrams using 
the FORM 3 |1(J| versio r@ of the package MINCER [TTJ- In some diagrams, e.g. with four 
external ^-fields where two external momenta are set to zero, this leads to IR divergences 
which mix with the UV ones in dimensional regularization. Another convenient method to 
compute renormahzation constants has been suggested in [12] and elaborated in the context 
of three-loop calculations in [J3]. The idea is an exact decomposition of all propagators using 
an auxiliary mass parameter M 2 : 

1 1 - p 2-2q-p-M 2 1_ 

(q+p) 2 q 2 -M 2 q 2 -M 2 (q+p) 2 ' ( ' 

3 For an n-point vertex V or a mass m the anomalous dimension is denned as = — /i 2 dl ^ v ' or 
7 m = — /1 2 dll ^i n ■ For a field we take the inverse renormahzation constant. 
4 The program can be downloaded from 

http : //www.nikhef .nl/~ f orm/maindir/packages/mincer/mincer .html 
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where p is a combination of external and q of internal momenta. This can be done recursively 
until the power in the denominator of the last term is high enough for this contribution to 
be finite, e.g. 

1 _ 1 -p 2 - 2q-p (-p 2 - 2q-p) 2 M 2 

(q+p) 2 ~q 2 - M 2 + {q 2 - M 2 ) 2+ (q 2 - M 2 f ~ {q 2 - M 2 ) 2 

M 2 (M 2 + 2p 2 +4q-p) (—p 2 — 2q-p — M 2 ) 3 1 ^ ' 

+ ( ?2 _ M 2)3 + ( g 2 _ M 2)3 (q+ P y 

As the result is independent of M 2 we can omit the contributions ~ M 2 in the above decom- 
position as long as we introduce counterterms into the Lagrangian to cancel M 2 -dependent 
subdivergences. In our case only a term ^-5Z^^ A'J 1 A afJ ' and ^-5Z^f^ $t<3> are possible!! 
The first one is not gauge invariant but this does not matter as it is only used for the can- 
cellation of subdivergences which works nevertheless. This method effectively amounts to 
introducing the same auxiliary mass parameter M 2 in every denominator of propagators and 
all possible M 2 -counterterms. We expand in the external moments!! and arrive at massive 
tadpole diagrams with one scale M. Due to the auxiliary mass no IR divergences can ap- 
pear while the UV counterterms which we are interested in (that is the ones without any 
dependence on the auxiliary mass M) will stay untouched. 

For the calculation of massive tadpoles we have used the FORM-based program MATAD [T3] . 
Where possible, i.e. for the propagators and three point functions, we have employed both the 
MINCER and the MATAD setups which served as an extra check. To generate the diagrams 
we used QGRAF [T5] and to compute the colour factors the FORM package COLOR [IB] ■ 

An important aspect of calculations such as these is the proper treatment of 75. As is well 
known, a naive treatment of 75 can be applied if it only appears in an external fermion line. 
In fermion loops we have to be more careful. In four dimensions we define 

75 = ztVtV = £j wVVW with e 0123 = 1 = -e 0123 . (20) 

In Fig. 1 to 5 we show a few diagrams that had to be calculated for the various ingredients 
of our final result. In order to have a contribution from a fermion loop with one 75 in it 
at least four free Lorentz indices or momenta on the external lines of the minimal subgraph 
containing this fermion loop are required. These can be indices from the gluon vertices or the 
internal momenta from other loops which act as external momenta to the minimal subgraph 
containing the fermion loop in question. External momenta of the whole diagram can be set 
to zero as the renormalization constants in the MS-scheme do not depend on those. 

Consider for example one of the fermion loops in FigJ2](c). The momenta on the two external 
$-legs can be set to zero. Then we have two indices from the gluon lines attached to our 
fermion loop and one loop momentum going through the two gluons and acting as an external 
momentum to the subgraph containing only our fermion loop. This is not enough to have a 
non-naive 75 contribution from this graph. 

For this reason diagrams like FigQ] (a,b,c,d), FigJ2] (a,b,c), Fig|3] (a,b), FigH] (a,b,c), FigJS] 
(a,b) can be treated naively. FigJ3] (c) has enough indices and momenta but no 75 in it. 
Diagrams like Fig{T] (c,e), FigJ5] (b), Fig|3] (b), FigH] (a,c), FigE] (a,c) are zero because of an 
odd number of 7-matrices in at least one fermion loop. And diagrams like Fig{T](b), FigJS] (b), 

5 Counterterms ~ M that would arise for fermions cannot appear because we have no M in the numerators 
of propagators. The ghost mass term ^SZffi c a c a does not appear because of the momentum dependence 
of the ghost-gluon- vertex. 

6 This method also works in massive theories. In this case we expand in the physical masses as well. 
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Fig|3](a), FigU](b), Figl3(b) are zero because of their colour structure. The only problematic 
type is FigfT] (f) which fortunately only contributes a ^ pole and can therefore be treated as 
described in [T7]. We use the fact that 75 anticommutes with every other 7- matrix in four 
dimensions and that 75 = 1. Then we apply relation (|2U|) for the case when one 75 remains 
on each fermion line. The two e^po can be rewritten as a combination of metric tensors 
which can be handled in dimensional rcgularization. The error we make with this treatment 
is of order e and does therefore not affect the pole part of our result. 






Figure 1: Some diagrams contributing to Z\ 
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Figure 2: Some diagrams contributing to 
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4 Results for the /3-functions 



First we give the results for the three-loop /3-functions of couplings A, y t and g s with the gen- 
eral gauge group factors for the strong interacting sector. Below C F and C A are the quadratic 
Casimir operators of the quark and the adjoint representation of the corresponding Lie alge- 
bra, d R is the dimension of the quark representation, T F is defined so that T F 5 ab = Tr (r a T b ^j 
is the trace of two group generators of the quark representation For QCD (colour gauge 



7 For an SU(iV) gauge group these are d R = N,C A = 2T F N and C F =T F (N-j r ). 
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Figure 4: Some diagrams contributing to Z\ 




group SU(3)) we have C F = 4/3 , C A = 3 , T F = 1/2 and d R = 3. Furthermore we denote the 
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(21) 



number of fermions by n f = Hi + 1. 
=12 A 2 + 2d R y 2 X - d R yf, 

/?f = _ 156 A 3 _ 2Mr y 2 x 2 _ y 4 A + g dBl/ 6 

+ 10CVd fl9 2 y 2 A-4C F ^ ffs V, 
=A 4 (3588 + 2016C 3 ) + 291d H y t 2 A 3 + y?\ 2 (™d R + 252C 3 d* - 3bd 2 

+ yfA y-^d a - 66( 3 d R + 804J + yf (yd fl - 12C 3 d fl - ^4 

(ont: \ 
— C F d R - 324( 3 C F d R ) 

+ f£v* t (-^C F d R + 60( 3 C F d R ^ + gffi\ (-^f-C F d R + 77C A C F d R 

-16n f T F C F d R + 72( 3 C 2 F d R - 36C 3 C A C F d R ) + g%f (^-C 2 F d R 

+48T F C F d R - ^-C A C F d R + 10n f T F C F d R - 48( 3 C 2 F d R + 24( 3 C A C F d R 

The purely A-dependent parts of eq. (|21|) have been known for a while |18p i9j. the full one- 
loop and two-loop result are in agreement with [221121] (for an SU(3) colour gauge group). 
It has also been a useful check for our setup to see that the same result can be derived from 
the four-<3?i vertex, the four-$2 vertex and the ($|$|$i$2) vertex. 

f =rf(! + irf»)-3<W, 

" \2 «„.2\ 1 „,4 9 



--3\ z -6yf\ + y?[---d R 



6C F + \c F d R \ + gt (~C 2 F - ^C A C F + ™n f T F C F , 



,a 2<2 /285 45 , \ 4 , /63 45 , 
18 A 3 + y 2 \ 2 - -d R j + y?\ + -d R 

345 . 9 „ 107 , 3. , 39 2 



+ V t {-32+Z<» + -n d «+2<» d « + l6 d *J (22) 

(57 81 
-C F + -C F d R 

+ gtvt (t^C 2 f - l -^-C 2 F d R + 2hT F C F + 7 -^rC A C F + 7 -]-C A C F d R 
V lb 8 lb 4 

33 27 

-— n f T F C F - An f T F C F d R - 27( 3 C 2 F + 18( 3 C 2 F d R - —Q 3 C A C F 

-% 3 C A C F d R ) + {~C^ + ^-C A C 2 F - ^C 2 A C F + 46n f T F C 2 F 



+^-n f C A T F C F + ^-n 2 T 2 C F - 48( 3 n f T F C 2 F + 48( 3 n f C A T F C F 



The one-loop and two-loop part of this result have been found before in [2D3E2] (for d R = 3, 
T F = i) and the contributions of order g 2 , gf, y 2 g 2 and gfy t to Z yt have been successfully 
checked against [23]. In this reference the calculation of Zm ± has been performed in the 
broken phase of the SM with a massive top quark. When comparing these two results one 
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has to take into account that in the broken SM the top quark mass is to be renormalized as 
a product yt ($2 + so * ne corresponding top quark mass renormalization constant is 

Zl , - 7 7^ - Z i 



7 (2t) „(2t) 
^2,L ^2,R 



Again the setup could be checked for consistency by using the renormalization of the four 
different vertices t-t-$2, t-t-<J>2> t-b-$i, t-b-$| for the calculation. 

ffl 2( , 2 



2T F g 2 y 2 + g\ (~C? A + 2n f T F C F + ^-n f C A T F , 



g 3 ' * sat ' 3S V 3 A ' ' ' 3 

= 4- n 2 7, 4 T-t„ 4 Zt„^ - n 4 7, 2 rsr„r7„ + i?r?.T^ ( 23 ) 



+ g^yt ^-T F + -T F d fl j - g*yf (3T F C F + 12C A T F ) 

V 108 lT 54^ n / A J F 

— ~9~ n/ F ~~ 27 / A ^ 

In the case of y t = this is in agreement with the well-known result |24tl25|. The one-loop 
and two- loop parts of eq. (f2"3"j) are known from [2"6Tf3"T] , the term oc gfy 2 can be found in [25] 
and the full three- loop result has been computed in [32] (for d R = 3, T F = ^). 
The /3-function describing the running of the "mass" parameter m 2 in eq. © can be computed 
from the renormalization constant of the local operator := 3>^<|>. An insertion of 02$ into 
a Green's function, e.g. with two external ^-fields, is renormalized as [O2*] = ^$2 02$ where 
[02$] is the corresponding finite operator. From [02$] = Z m iO\^§ c and O^ 6 = Z < 2^C>2<$> it 
follows that 

Z m * = (Z^y 1 (24) 
This yields the following contributions to j3 m 2 : 

PI 

=6X + d R y 2 , 

B {2) Q 

— = - 30 A 2 - 12d R y 2 X - -d R yf + 5C F d R g 2 y 2 , 

/3 (3 2 } , 99 9 9 

m z 2 V 4 

611 

Te 

+ 5s \ 2 A (-153C F d R + 144C 3 F d fl ) 
+ 9 s V(^C F ^-90C 3 O F d fl 



99 / 333 

=1026 A 3 + — d R y 2 X 2 + yfX ( —d, 

+ ^(-^-d R + 24dl + 15C 3 d R ) (25) 



/ iin 77 \ 

<? s \ 2 y—C 2 F d R + —C A C F d R - 8n f T F C F d R + 36( 3 C 2 F d R - 18( 3 C A C F d R ) 



The one-loop and two-loop parts of this result are in agreement with [2U] where they have 
been computed before. The purely A-dependent part can be found in JT8JQ35]. For d R = 3 
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and T F = \ (QCD) we get the following results: 

=l2\ 2 + 6y 2 \-3yf, 
/3i 2) = - 156 A 3 - 72y t 2 A 2 - ^ yfX + 15y t 6 + 40 5 2 y t 2 A - 16 5 2 y t 4 , 



=A 4 (3588 + 2016Cs) + 873 y 2 t A 3 + y t 4 A 2 + 756C 3 ) 

+ yfX - 198C 3 ) - yf + 36C 3 ) + 5 2 y 2 A 2 (-1224 + 11526) 

+ g 2 3 yt\ (895 - 1296C 3 ) + 9^ (-38 + 240C 3 ) 
+ ^ 2 A ( ^ - 32 % - 48C 3 ) + atvt (~ + 20n f + 32C 3 



(26) 



V 3 1 " J " si " V 3 

To get an idea of the size of these contributions and therefore the significance of our calculation 
we evaluate f3 x at the scale [i = M z (with an assumed Higgs mass of 125 GeV and n f = 6) 
which yields a value of /3 A ~ (—0.01) at one-loop level. The two and three-loop contributions 
change this result by ~ 1% and ~ (—0.04)% respectively. To estimate the importance of the 
individual terms we introduce the labels 

p — 9s y _ Vt j _ A 

g a (fi = M z )' y t (fi = M z y A(/i = M z ) y > 

and get 

1 loop / \ 1 loop 1 loop 

-8.5 G 2 Y* +5.0 Y 6 +3.1 G 2 LY 2 \ 

2 loop 2 loop 2 loop 



(28) 



+ I 7.9 G 2 Y 6 -4.8 Y 8 -5.3 L 2 Y 2 -3.1 G 2 LY A 

3 loop 3 loop 2 loop 3 loop 

- 2 . 5 G A Y 4 + 2 . 6 G 4 2 - 1 . 7 L 3 | 10~ 5 

3 loop 3 loop 2 loop 

+ I -7.5 < LY A +7.8 L 2 Y 4 -6.6 LY* +L1 G 2 L 2 : 

\ 2 loop 3 loop 3 loop 3 loop 

+ I 5.7L 3 y 2 +5 : 9L 4 I 10~ 7 . 

\ 3 loop 3 loop / 

We see that the decrease of the effective four-<3? coupling with increasing energy is induced 
by top quark loops. Without quarks there would be an increase. It is also worth noting that 
the individual contributions at three-loop level are much larger than the overall effect due 
to huge cancellations. Consider for example the five numerically largest three-loop terms at 
H = M z : 

r 7.9 G 2 Y 6 - 4.8 Y 8 - 3.1 G 2 LY i - 2.5 G 4 Y 4 + 2.6 G*LY 2 ) 10" 5 . 



The total contribution from these terms is by almost two orders of magnitude smaller than 
the size of the largest one. 
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For the top- Yukawa /3-functions we find 



4" 9 2 „ 2 

^S- -3 A 2 - erfA - 6rf + 18<« + g', (-?f + | % 

+ 8«A-f 9 ,V + ^(^-f n; -114 Ci 



+ 



fi / , 2216 140 2 160 A \ 
«? ("1249 + ^ + —^ + — Can,) 



As has already been mentioned above the A-corrections are negligible here. Evaluating /? 
at the scale \x = Mz (with an assumed Higgs mass of 125 GeV and n } = 6) we get a value 
of ~ (—0.023) at one-loop level which means a decrease of y t and therefore the top mass 
with increasing energy. This is due to the QCD corrections. In the absence of QCD the 
opposite would be the case as we can see from the term oc y 2 in The two and three-loop 
corrections are ~ 16.6% and ~ 0.7% with respect to the one-loop result and so quite high 
compared e.g. to the case of f3 gs discussed below. We use again the labels (|2"7|) to get an 
impression of the individual terms: 



Py t \n=M z = \ r 3.7YG 2 +1.3Y^ W~ 2 + -4.7JG 4 , | 10" 

1 loop 1 loop / \ 2 loop 

+ [ 9.8Y 3 G 2 -2.5YG e -2.1Y 5 ) 10~ 4 

2 loop 3 loop 2 loop 

9.3 Y 3 G 4 -3.1 LY 3 -2.5 Y 5 G 2 | 10~ 5 (30) 

3 loop 2 loop 3 loop 

+ | 6.0y 7 +3.0Ly 5 +2.3L 2 y 1 10" 6 

3 loop 3 loop 2 loop 

I 3.9L Y 3 G 2 J 10~ 7 + I -1.2^L 3 y I 10~ 8 + I 8.5 L 2 Y 3 

3 loop / \ 3 loop / \ 3 loop 



For the strong coupling we get 



2/11 1 

— =9 h -n f 1 

9s 93 V 2 T 3 1 



9 2 a V 2 t +9t(-51+ 1 4n f ), 



9s V 3 

M 3) 15 2 4 42 R / 2857 5033 325 2 

i-2!_ =— oV - 20 g 4 y 2 + o 6 + n f n 2 

g s 2 ysUt ysyt ys V 4 36 7 108 / 

In order to numerically compare the higher order corrections to the above /3-functions we also 
give the evaluation of (3 ga at the scale [i = Mz and with n } = 6. The one- loop contribution 
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is ~ (—0.04) to which the two and three-loop calculations give corrections of ~ 3.7% and 
~ (—0.02)% respectively. Very small A-corrections to /3 9s do not appear until four loops. 
With eq. (|2T|) and the above assumptions we get 



AJm=m z = | r±2£l) 10-2 + d^£t 1 10-3 

1 loop / \ 2 loop 

-6.9 Y 2 G 3 +1.7 G 7 -1.3 Y 2 G 5 I 10~ 5 + I 3.1 Y 4 G 3 | 

2 loop 3 loop 3 loop / \ 3 loop 



(32) 



The running of the m 2 parameter is given by 

/? (1 2 } 

\ =6A + 3y t 2 , 



,(2) 

30X 2 -36y 2 \-^y? + 20^ 
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«« \2 o«„,2\ z ' „,4 , on.2„,2 



^=10^ + ^ 2,2+4^^1 N (33) 



^-1026 A a + =^ ^ + yf A + 216( 3 

+ V & t + 45 C3 ) + 5s \ 2 A (-612 + 576C 3 ) 

+ - 360C 3 ) + 9% 2 - 16n, - 24 C3 ) . 

Again we evaluate /3 m 2 at the scale \i = Mz (with an assumed Higgs mass of 125 GeV and 
n f = 6) and get a value of ~ (+0.023) at one-loop level which means an increase of m 2 at 
higher energy scales. The two and three-loop corrections are ~ 2.9% and ~ 0.32%. With the 
labels (|27|) we can estimate the contributions of the individual terms: 

\^= Mz = 10~ 2 + I 5^+1. lG 2 ^ 2 I 10~ 3 + I -2.4y 4 -1.9Ly 2 I 10~ 4 

\ 1 loop / \ 1 loop 2 loop / \ 2 loop 2 loop / 

+ I 9.4 G 4 Y 2 -7.0 G 2 Y 4 +3.3 Y 6 -2.4 L 2 +1.1 LY A I 10~ 5 (34) 

\ 3 loop 3 l 00 p 3 loop 2 loop 3 loop / 

1.0 G ' 2 LY' 2 J 10~ 6 + I 7_1L^ +7.0 L 2 Y 2 | 10~ 7 

3 loop / \ 3 loop 3 loop 



5 Results for the anomalous dimensions 

In this section we give the anomalous dimensions of the physical fields in this setup. Note 
that because of the SU(2) symmetry $1 and $2 must have the same anomalous dimension 
72 . The same holds for the left-handed part of top and bottom quarks: I2 L = ^2 L- ^or the 
quark flavours q which do not participate in the Yukawa interaction there is no difference 
between the left- and right-handed part as they are renormalized by the same Z-factor: 
T2 = 72 l = ^2 w This also applies to the right-handed part of the bottom quark: 72 R = 7|. 
All these relations have been tested explicitly during our calculation which provides a nice 
additional check. 
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7 « (1) = 5 2 C F (l-£), 

72 (2) =9s (-\c 2 F + ^-C A C F - In fTpC p - ^C A C F + -^ 2 C A C F , 

7 ^ (3) =6T F C F gfy 2 + gt (§C? - + ^g^C F + 3n f T F C F (35 ) 

-^n.C^CV + ™n 2 T 2 C F + 12C 3 C A C F - yCaC^ - ^C 2 A C F 
+^n f C A T F C F - \^C 2 A C F + ^ 2 C 2 A C F + ^ 2 ( 3 C 2 A C F - ^ 3 C 2 A C F 

For y t = this is in agreement with the well-known QCD result [25J. The renormalization 
constants for t,b and q can also be found in [23] up to order g® and g\ y 2 where the calculation 
has been performed in the broken phase of the SM. 

9(1) , 1 2 

l2,L =72 + g ^. 

72? =7l (2) - yt 4 Q + ^)-2c,^; 

72?? =7^ (3) " f l/ t 2 A 2 + + (-| + f d fl - -d R + - C 3 ) (36) 

/ 1 Q 5 

+ (- y C 2 F + yC A C F - ^n / T F C F + 6( 3 C 2 F - 15( 3 C A C F 



2„.2 

t 1 



72,R —72 + Vt 1 



l2 { S=l2 i2) -yt(\ + ld R yAC F9 2 y 2 , 



t(3) ,(3) 33 2 2 4 6 / 33 53 3 ,j2 1 qa \ 

7 2 ,V =72 2 Vt + 12 Vt A + Vt 16 + ~8~ * " 4 + 3< ^ 3 J (37) 

+ 5s V (W + ^C F d fl + l2( 3 C F d R 

+ 9tv 2 i (-^-Cj + ^-C A C F - 3n f T F C F + 12( 3 C 2 F - 30( 3 C A C F 



-A 7 ,2 

72 — a R tit j 

7 * ( 2 ) =6 A 2 - Hd B y f + 5C F d fl g 2 yl 

7 *(3) = _ 36A 3_^ y 2 A 2 + 1Mfly 4 A 

+ 2/ t 6 + 64 + 3Cs^) + g 2 yf {^C F d R - lK 3 C F d R ^ 

f 119 77 \ 

+ 5^ t 2 f -— C 2 F d R + - 8n f T F C F d R + 36C 3 C 2 d« - 18( 3 C A C F d R ) , 

(38) 



13 



The purely A-dependent part of this has been computed before in |18yi9|. 
7 2 9(1) =i{-\c A + \n f T F -\iC^ 

7 | (2) = - 4T F + gt {- 2 ^Cl + AufTpCp + 5n f C A T F - + 1^ ) , 

72 (3) =9 2 sVt (9T F + 7T F d R ) - g% 2 (gT f C f + ^C A T f 

JA n ^C A T 2 F + hzCl + 24( 3 n f C A T F C F - 18( 3 n f C 2 A T F - ^C 3 A 
9 2 lb 

This is also in agreement with [25] for y t = and with [23] for y t = 0, £ = 0. For d fl 
T F = \ (QCD) these results £1X6 clS follows: 

7| (1> =| (1 " rf, 

.4rf, ? + I ( + ^ - ^n, + |„ ? - 26f 3 - ^ { 



9(1) , 1 2 

7 2 ,l =7 2 + g ft. 

t(2) 9 (2) 5 4 8 2 2 

7 2 ,l =7 2 - g 2/4 ~ 3 9aVt ' 

t(3) 9 (3) 33 2 2 . c 4\ , 6 , 3 a 

7 2 ,l =7 2 - -J Vt A + 6 ?/t A + Vt ( 6 + 2& 

+ 3 S V f y + 24Ca) + S 4 y t 2 ( y - % - — C 3 



i2,R —'2 + y< ' 

t(2) ,(2) 19 4 16 2 2 

7 2 ,V =7 2 4 y * ~Y 9 ' Vt ' 

t(3) ,(3) 33 2 2 4 6 /177 

7 2 ,V =7 2 - yy t A + 12 y t A + y° ^— + 3C 3 

2 4 /35 \ 4 2 /25 „ 296 



+ g 2 y? ( - + 48C 3 ) + gtv 2 ( y - 2% - — Cs ) , 



7 2 * (1) =3 yl 

72 (2) =6A 2 -f y 4 + 20 5s \ 2 , 

7* < 3 > = - 36 A 3 - ^ y 2 A 2 + 45 y?\ + + 9 C 3 

+ ff^ 4 (y " 72 C3 ) + ^ (jf - 16% - 24 C 3 
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i r 2 4 yl 4 2 



(44) 




3429 
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6 The evolution of the quartic Higgs coupling 

The quartic Higgs coupling A is of special interest as it is directly related to the Higgs mass 
Mh- If we assume that the SM is valid up to some high energy scale A, then the value of 
Mh should meet the constraints 



Here the upper limit is related to the well-known fact that the running Higgs self-coupling 
develops a Landau pol^l if Mh is large [3"3Tf3"5] . For A = Mpi anc k = 10 18 GeV the estimated 
value of m max is around 175 GeV [55H55] , which is already excluded by experiments carried 
out at the LHC and the Tevatron. 

The lower limit m m i n follows from the requirement of the vacuum stability |37H39| . In order 
to find m m i n one should construct the effective Higgs potential V[4>] including radiative cor- 
rections and sum possible large logarithms using the standard method of the Renormalization 
Group (for a review see, e.g. [10]). Once this has been done, the condition that the potential 
V[4>] does not develop a deeper minimum in addition to the standard one for all values of 
4> < A fixes m min . 

In our analysis we will use a simplified approach for finding m m j n , namely the requirement 
that the running coupling constant A(/x) stay non-negative for all \i less than A. It has been 
shown in |34U41j that the simplified approach is essentially equivalent to the one based on the 
use of the effective potential provided the instabilty of V[(f>] can only happen at <j) 3> Mz- 

(3) 

In this section we investigate the effect of the three-loop result /3 X on the running of A and 
therefore its effect on the stability of the electroweak vacuum in the SM. For this we also 
include the electroweak contributions up to the two-loop level. The two-loop /3-functions for 
the SM gauge couplings have been derived in [28H3T] . The two-loop results for the Yukawa- 
couplings and A can be found in [20H22] . Now we add the three-loop results derived in the 
previous section and investigate the effect this has on the evolution of our couplings (for a 
recent similar analysis, using the two-loop running, see, e.g. [9"ll42tl4"3"] ) . 

To find starting values for the running of the couplings we should account for the fact that 
the physical parameters (e.g. pole masses) are related to the ones in the MS-scheme in a 
non-trivial way (see e.g. [44H46| ) . For example eq. (|15|) is only valid at tree level. For the 
higher order corrections we take the electroweak ones at one-loop and the QCD ones at two- 
loop level from ( J 1. 15 . These matching relations depend on the exact values of a a {Mz) and 

8 This is true in the one-loop approximation. At two loops the Landau pole is replaced by an ultraviolet 
metastable fixed point with the resulting fixed point value of A being outside the weak coupling region. 



< M H < m, 



max • 
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the pole mass Mt of the top quark and of course the mass of the Higgs boson Mjj. For the 
latter we consider the cases Mjj = 124 GeV and Mh = 126 GeV. For the other two we use 
the values 

a s {M z ) = 0.1184 ±0.0007 and M t = 172.9 ± 0.6 ± 0.9 GeV g7j. (45) 

One should also keep in mind that the matching relations themselves receive contributions 
from not yet known higher order corrections . The corresponding uncertainty in m m i n has 
been estimated in [UH2] and found to be about 2 GeV. 

Fig. [U] shows the evolution of A in this framework up to the Planck scale. To estimate the 
dependence of the A-running on the parameters a a {Mz) and Mt we give the shifted curves 
for X(fi) when we change these parameters by ±<r as given in eq. (|43|) l^l As the two- and 
three-loop curves are very close together we zoom in on the region where A crosses over to 
negative values in Fig. [7J In this plot we give the a s -uncertainty for the two- and three-loop 
curves to compare between this uncertainty and the shift from two to three loops. Note that 
there is a considerable difference between Mh = 124 GeV and Mjj = 126 GeV which means 
that the evolution of A is very sensitive to the value of the Higgs mass. Given a fixed value 
for Mh the largest uncertainty lies in the exact value of the top mass. The second largest 
uncertainty comes from a a . The total effect due to the three- loop part of the /3-functions is 
somewhat smaller than latter as can be seen best in Fig. [71 Still, it is worthy of note that 
the three-loop corrections to the /3-functions presented here enhance the stability of the SM 
electroweak vacuum. 

The smallness of the three-loop correction to /3 A seems to be somewhat coincidental as the 
aforementioned cancellations of individual terms in /3® depend strongly on the value of 
Mh- Finding a Higgs with a mass of 124 to 126 GeV would therefore mean an excellent 
convergence of the perturbation series for /3 A . Another intriguing consequence of a Higgs 
mass in that region is the uncertainty whether A becomes indeed negative at high scales 
or not. If we take e.g. Mh = 126 GeV and a s = 0.1184 and decreasqij the top mass from 
M t = 172.9 GeV to M t = 171.25 GeV (M t = 171.16 GeV without the three-loop corrections), 
then A stays positive up to the Planck scale Mpi anc k in our framework. The same effect can be 
achieved for M H = 126 GeV and M t = 172.9 GeV by increasing a s = 0.1184 by 6.5 a as (7 a as 
without the three-loop corretions). A combined scenario for Mh = 126 GeV would be a shift 
of M t = 172.9 GeV by -la Mt = -1.5 GeV and of a s = 0.1184 by +la as = +0.0007 GeV 
which would also make A positive up to the Planck scale. 

Thus, we conclude that at present no definite answer can be given to the question whether 
the SM vacuum is stable all the way up to the Planck scale or not. If indeed a SM Higgs 
boson is found with a mass of 124 to 126 GeV, this is a good motivation for determining 
a s (Mz) and Mt as accurately as possible as well as calculating the SM /3-functions to the 
highest achievable accuracy. 

7 Conclusions 

We have computed the three-loop corrections to the evolution the top- Yukawa coupling, the 
strong coupling and the quartic Higgs self-coupling in the unbroken SM with the numerically 
small gauge coupling constants g± and g<i and all Yukawa couplings except for y t set to zero. 

9 In order not to make the plot too crowded these shifted curves are only given for the two- loop result. 
The difference to the three-loop result is similar to the one between the two- and three-loop curves for 
a„(M z ) = 0.1184 and M t = 172.9 GeV. 

10 Smaller values for Mt or larger values for a B increase the stability of the vacuum. 
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Figure 6: Evolution of A with the scale /i: 2 loop (dashed, blue) and 3 loop (continuous, red) 
results; Uncertainties with respect to the two-loop result: ±lcr aa , ±lo"M t (dotted) 

The implications of our calculation on the stability of the electroweak vacuum in the SM can 
be summarized as follows: 

• The total effect of the three-loop terms is relatively small which is not self-evident as 
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Figure 7: Evolution of A with the scale /i: 2 loop (dashed, blue) and 3 loop (continuous, red) 
results; Uncertainties with respect to the two- and three-loop results: ±l<7 as (dotted) 



the individual terms in f3 x are much larger than the final value due to significant 
cancellations for a Higgs mass in the vicinity of 125 GeV. 
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• The evolution of A is very sensitive to the values of the Higgs mass, the top mass and 
a 3 (Mz)- If we take e.g. Mh = 126 GeV and decrease the top mass by about 1.7 GeV, 
then A stays positive up to the Planck scale (a similar observation has been made 
in 02]). This is a very good motivation for high precision measurements of a s (Mz) 
and Mf. With the latter values known more precisely, the account of the the three- 
loop effects in the evolution of the quartic Higgs self-coupling would be essential in 
considering the problem of the stability of the electroweak vacuum in the SM. 

• In this context it may also be useful to calculate both the electroweak contributions to 

and at three-loop level as well as the matching of experimentally measurable 
on-shell parameters and MS-parameters to a higher accuracy. 

We thank Luminita Mihaila, Jens Salomon and Matthias Steinhauser for useful discussions 
and informing us on the results of [32] before their publication. We thank Johann Kuhn for 
valuable comments and support. Last but not least, we are grateful to Fedor Bezrukov and 
Mikhail Shaposhnikov for gently drawing our attention to the subject, numerous discussions 
and providing us with their version of a Mathematica package to perform one-loop matching 
in the SM. 

In conclusion we want to mention that all our calculations have been performed on a SGI 
ALTIX 24-node IB-interconnected cluster of 8-cores Xeon computers using the thread-based 
|48] version of FORM [TU] . The Feynman diagrams have been drawn with the Latex package 
Axodraw [39] . 

This work has been supported by the Deutsche Forschungsgemeinschaft in the Sonder- 
forschungsbereich/Transregio SFB/TR-9 "Computational Particle Physics". 

Note added. More elaborated analyses of the vacuum stability which take into account the 
three-loop running as obtained here and two-loop matching corrections can be found in recent 
works [SDH53]. 
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